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We examine the dynamical magnetic response of the underdoped cuprates by employing a phe- 
nomenological theory of a doped resonant valence bond state where the Fermi surface is truncated 
into four pockets. This theory predicts a resonant spin response which with increasing energy (0 
to lOOmeV) appears as an hourglass. The very low energy spin response is found at (tt,tt db S) and 
(tt ± 5, 7r) and is determined by scattering from the pockets' frontside to the tips of opposite pockets 
where a van Hove singularity resides. At energies beyond 100 meV, strong scattering is seen from 
(7r, 0) to (71*, 7r). This theory thus provides a semi-quantitative description of the spin response seen 
in both INS and RIXS experiments at all relevant energy scales. 
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Introduction: Neutron scattering studies of the mag- 
netic properties of underdoped cuprate superconductors 
have revealed an unusual hourglass pattern in the spin 
excitation spectrum that persists into the normal state 
PP. This spectrum which is centered on (71", 7r), can be 
divided into three energy regions. At low energies the 
weight is shifted to nearby incommensurate wavevectors, 
peaking along the crystal axes. With increasing energy 
the weight is more uniformly distributed around (-zr, tt) 
and contracts to form a resonance centered on (71", tt). At 
still higher energies a uniform ring appears evolving away 
from (tt, 7r). Recent RIXS experiments [14 have explored 
this high energy region further. 

A phenomenological theory for the underdoped pseu- 
dogap phase by Yang, Rice and Zhang (YRZ) [2] has 
had considerable success in reproducing many electronic 
properties [3]. In this letter we examine the spin spec- 
trum within this theory and show that key features of 
the experiments at all three energies are reproduced. We 
begin with a derivation of the YRZ ansatz starting from 
a t- J model rather than from the overdoped Fermi liq- 
uid in the original paper, following a recent suggestion 
by P. A. Lee [4 which in turn is based on a small mod- 
ification of Ref. [5]. An RPA form is used for the spin 
response similar to that employed by Brinckman and Lee 
in their study of the spin resonance in the superconduct- 
ing state of overdoped cuprates. The itinerant descrip- 
tion in our approach differs from the often used scenario 
of strong and slow stripe fluctuations which would show 
up as incommensurate quasi-elastic peaks in the mag- 
netic response [7]. We compare our results with data 
sets gathered using both x-rays and neutrons. 

YRZ Spin Response: The YRZ ansatz, as originally 
conceived, was an ansatz for the single particle Green's 
function (GF) of the underdoped cuprates. The Fermi 
surface associated with it is truncated and composed of 
four nodal pockets (see Fig. 1). The area of the pock- 
ets is proportional to the level of doping, x. This GF is 
also characterized by lines of Luttinger zeros which coin- 




FIG. 1: a) The Fermi surface for hole doping, x = 0.12. 
Hole pockets are marked in red (solid) while the lines of 
Luttinger zeros are blue (dashed). Also marked is a nesting 
vector Q = (0.5,0.375) (in reciprocal lattice units), connect- 
ing the tip of a pocket to the frontside of another pocket. 
Inset: A plot demonstrating that at the tip of the pocket 
there is a saddlepoint in the superconducting quasi-particle 
dispersion and hence a van Hove singularity. Contours are la- 
belled in meV. The parameters used here are t(x) = 70meV, 
t'( x ) = -0.18t(x), t"(x) = 0.12t(x), A = 0.34t(x) and 
Asc = 0.05£(x). b) The real parts of Xyrz(Q,oj) and 
Xo(Q,cj) vs. uj. c) The imaginary parts. 



cide with the magnetic Brillouin zone (BZ) or Umklapp 
surface [6] (see Fig. 1). The ansatz was inspired by an 
analysis of a system of weakly coupled Hubbard ladders 
where a similar phenomenology was found to hold [8J. 

To extend the YRZ ansatz from the single particle GF 
to the spin response, we first elucidate the connection 
between YRZ and the slave boson (SB) treatment of the 
t-J Hamiltonian. SBs provide a natural RPA-like form 
to the spin response and we intend to adapt this to the 
assumptions of YRZ. For this purpose we then write the 
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t — J Hamiltonian as 
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The Hamiltonian is divided into terms involving near- 
est neighbour (NN) hopping, i^ £ nn , terms involving next 
nearest neighbour (NNN) hopping (and beyond), H™ nn , 
and a spin-spin interaction, Hj H . We choose this sepa- 
ration because of the focus the YRZ ansatz places upon 
the Luttinger zeros found at the magnetic BZ. The near- 
est neighbour dispersion, £o(&) = — 2t(x)(cos k x + cosk y ) 
(for the definition of t(x) see [2j[T0]) vanishes on this line 
while that of the NNN hopping does not. (We show in the 
supplementary material [10] that one can derive a YRZ- 
like propagator while treating NN and NNN hopping on 
the same footing.) We now subject H^ n + Hj H to the 
standard slave boson mean field treatment (leaving H™ nn 
to later). We thus factor the fermions, cj^, into spinons, 



fj a and holons, bi via c\ a = fj a bi, where the spinons and 
holons are subject to the constraint ^ a fi a fia + b\bi = 1. 
At this level the spinon Green's function can be shown 
to be [9] 
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where = \A R (k)\ 2 /(uj + f (k)) and A R (k) = 

A (x)(cos k x — cos k y ). Here t(x) and A (x) are doping 
dependent parameters. The single particle Green's func- 
tion, G£, is given directly in terms of the spinon Green's 
function because we assume the bosons are nearly con- 
densed and so replace the boson propagator (b\(r)bj(0)) 
by gt(x): G%{uj, k) = g t (x)G{ r (u, k) (in slave boson mean 
field gt(x) = x [9 j; in the Gutzwiller approximation 
g t (x) = 2x/(l + x) 0). This is close to the YRZ form 
but differs in that the full dispersion in the denominator 
is replaced by the dispersion due to NN hopping. 

To bridge the gap between the SBMFT and YRZ, we 
then turn to the so far neglected NNN hopping, H^ nn . 
Treating this term in mean field theory (MFT) moves 
the Luttinger zeros off the magnetic Brillouin zone and 
so we instead use an RPA like approximation (see Fig. 
2) leading to 
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Here £'(fc) = — £t f {x) cos k x cos k y — 2t"(x)(cos2k x + 
cos2k y ) is the dispersion due to the NNN terms. The 
spinon propagator in this form now immediately gives 
the YRZ ansatz. An important consequence of the non- 
MFT treatment of the H™ nn is that spinons and holons 
are bound together. This binding distinguishes the YRZ 
ansatz from the standard mean field SB approximation 
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FIG. 2: RPA form of the YRZ spinon propagator in terms of 
SB propagators. 



which leads to an expanded Hilbert space with indepen- 
dent spinons and holons. A second consequence is the 
absence of an anomalous spinon propagator (or at least 
the coherent part thereof), consistent with an underly- 
ing RVB assumption that spin correlations are short- not 
long-ranged in the YRZ ansatz. 

This form (Eqn. [8| applies in the normal phase and 
can be generalized to the d-wave superconducting state, 
e.g. see [3]. Note that YRZ gives a two-gap description 
of the pseudogap phase with a separate RVB (Ao) and 
pairing A^c gaps. 

We now turn to the YRZ spin response. In slave 
bosons, neglecting the effects of NNN hopping, the spin 
response naturally takes on an RPA-like form [9]: 
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Here xo(^?k) is the bare particle-hole bubble for the 
spinons (including anomalous contributions) and J(k) = 
J(cos k x + cos ky). 

How now does our non-mean field treatment of H™ nn 
alter this? Its effect is two-fold. Firstly we no longer 
include a contribution to xo from the anomalous spinon 
Green's functions. And to determine how t nnn dresses 
the normal spinon Green's functions, we employ the same 
approximation that led to the YRZ ansatz. Namely we 
only allow diagrams involving vertices where the boson 
lines of the vertex are tied together. Under such a restric- 
tion, t nnn only dresses the individual spinon propagators 
making up the particle-hole bubble entering xo- I n this 
fashion the YRZ spin ansatz takes the form 
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where Xyrz is simply a particle-hole bubble made up of 
YRZ quasi-particles. 

In computing the spin response we treat J as a fit- 
ting parameter for each doping (and different from «/#). 
We do not expect the underlying mean field treatment 
to accurately treat the renormalization of J which is in- 
evitably doping dependent. In particular in the presence 
of strong scattering connecting the magnetic Brillouin 
zone boundaries, we expect J to be strongly modified. 
This is not merely a feature of the YRZ theory but is 
generic to slave boson flavoured theories. In Ref. 0, J 
had to be sharply reduced in order to produce an AF 
ordering transition at approximately the correct doping. 
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Results: We begin with the lower energy (uj < 
lOOmeV) spin response in the underdoped cuprates which 
has a universal hour glass shape [Tj [TTHT3] with strong in- 
commensurate response at low energies (i.e. uj ~ 2A$c) 
concentrated at four points, (7r,7r ± 5) and (jr ± 5, 7r). 
As energy is initially increased this response evolves in- 
wards towards (7r,7r) and simultaneously becomes more 
isotropic in its distribution about (7r,7r). Whether this 
inward dispersion reaches (7r, tt) is a function of the par- 
ticular cuprate being examined. As the energy is further 
increased this inward evolution is reversed and the re- 
sponse moves outwards from (7r,7r). With this outward 
dispersion, the response is more isotropically distributed 
about (7r,7r). 
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FIG. 3: Constant energy slices of the spin response for x = 
0.12 in the SC phase - the parameters used are the same as 
listed in Fig. 1, with J — 0.14meV for our theory. 

We see these general features in the constant energy 
scans of the q-dependent spin response presented in Fig. 
3 for the superconducting case. In this figure we have 
chosen parameters appropriate for the description of un- 
derdoped La2- x Sr x Cu04. We see at very low energies 
(0.05 J) the primary response is at (7r, 7r± S) and (7r±5, tt) 
with S = 0.167T. As the energy is increased there is a 
slight inward dispersion (5 decreases slightly) and the 



spin response is found circularly distributed about (7r, tt). 
This dispersion reverses at uj ~ 0.2 J and begins to move 
outwards. In this energy range the greatest response is 
found about (tt ± 5' , tt ± 5'). 

The response found at (tt, tt ± S) and (tt ± <5, tt) at 
0.05 JmeV can be directly ascribed to transitions between 
the fronts of the pockets and the tips of opposite pockets 
(see the vector Q in Fig. 1). In general the presence 
of the pockets in the YRZ theory allows for low energy 
scattering in a larger portion of the Brillouin zone than in 
theories where the spinon Fermi surface consists of four 
points coinciding with nodes of the SC order parameter 
(see Figs, lb and lc for a comparison of Xyrz and xo; 
Xo is the bare particle-hole bubble for the standard slave 
boson description of the spin response [9]). Moreover in 
the presence of a SC gap, the tips of the pockets see a 
saddle point in dispersion with a corresponding van Hove 
singularity thus further enhancing the low energy scat- 
tering. 
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FIG. 4: a) Hourglass dispersion of the resonance near 7r,7r. 
The thick black line is the position of the maximum in- 
tensity peak after integrating the numerical data over a 
strip of width 27r/25 along the parallel direction, averaged 
over sections of length 27r/33. Experimental data points 
(appropriately rescaled) are taken from [TJ [T5] : 'up' tri- 
angles, Lai. 90 Sr .ioCu04 [II] ; circles, Lai.87sBao.i25Cu04 
[16] : 'down' triangles, Lai.84Sro.i6Cu04 [13] : squares, 
YBa2Cu306.5 [E] and diamonds, YBa2Cus06.6 hayden. b) 
k-integrated spin response with and without a superconduct- 
ing gap. 

With increasing energy, the k-points with maximal in- 
tensity move inward towards (7r,7r), albeit in an uneven 
fashion (there is a sudden movement inward at 0.125 J) 
while at the same time becoming more isotropically dis- 
tributed about (7r,7r). This behavior is shared not only 
by underdoped La2- x Sr x Cu04 [11 but also its optimally 
doped counterpart [T2l [13] . It is also seen in stripe sta- 
bilized La 2 - x Ba x Cu0 4 [E] and YBCO [U US]. At an 
energy between 0.4J and 0.5J, the point of maximal in- 
tensity begins to drift outward from (tt, 7r), again a uni- 
versal feature of the magnetic response in the cuprates. 
We explicitly plot in Fig. 4a the evolution of the k-point 
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of maximal intensity as a function of energy, comparing 
its evolution with a number of cuprates. 

In the normal state, low energy spectral weight is found 
not just in directions parallel to the crystal axes but in 
the nodal directions as well (see Fig. 5). While parallel 
scattering still dominates at low energies, the response is 
less concentrated in such areas and weight does appear 
along the nodal directions (at least in the LSCO family) 

[mini. 



The paramagnon results from a particle-hole excitation 
from the lower band, cj_(k), to the upper band, cj+(k). 
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FIG. 5: x = 0.12 and uj — 7meV constant energy slices for 
the SC phase (left) and the normal phase (right). 

Underlying our calculations of the magnetic response 
is the assumption that itinerant quasi-particles (even if 
heavily dressed) can explain this response in the cuprates. 
While there is evidence that at least part of the spin re- 
sponse must be ascribed to localized spins |16j [19] , there 
is also evidence that impurities introduce local spins, e.g. 
Zn doped into YBCO [20] and earlier studies. The full 
cuprate magnetic response requires a mixture of the two. 
However one experimental feature of the spin response 
that points to itinerant quasi-particles is the depression 
of the k-integrated spin response at uj < 2Asc upon de- 
creasing T < T c . This behavior is seen in both the LSCO 
[TTHT3] and YBCO families and we see it as well in our 
calculations (Fig. 4b). We also see in Fig. 4b that our 
calculated integrated intensity has a two peak structure, 
with one peak at energies close to 0.05 J and one at en- 
ergies at « 0.12 J. This doubling of peaks is seen in near 
optimally doped LSCO [HEL3]. In underdoped LSCO at 
least the lower energy peak has been observed [TT] . 

Turning to high energies, uj > lOOmeV, we find the 
YRZ spin response is also able to explain key features 
in the spin response recently measured by RIXS exper- 
iments. In Fig. 6 we plot the spin response for ener- 
gies lOOmeV < uj < 300meV for two cuts in the Bril- 
louin zone. We see two features emanating from (0,0). 
One disperses towards (-zr, 0) as energy is increased (cor- 
responding well with the reported paramagnon-like exci- 
tation in the RIXS data of [14J on a variety of cuprates). 
The other, with a considerably greater spin velocity, 
evolves towards (71,71). This dispersing paramagnon ex- 
citation naturally appears from a two-band factorization 
of YRZ where the propagator can be written in the form 
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FIG. 6: The spin response for energies from lOOmeV to 
300meV, for cuts from (tt,tt) to (0,0) to (tt, 0) in the Bril- 
louin zone (same choice of parameters as previously). Also 
plotted are data points from [14]: circles, Ndi^Bai.sCusOe; 
squares, YBa2Cus07; diamonds, Ndi^Bai.sCusO?; 'up' tri- 
angles, YBa2Cu40s; 'down' triangles, YBa2Cus06.6- 

In conclusion we have shown that calculations of the 
magnetic response based upon itinerant YRZ quasi- 
particles reproduces key features reported in experiments 
on the spin response of the underdoped cuprates at both 
low and high energies in a satisfactory way. 
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Supplementary material 

In this appendix we give an alternate derivation of 
YRZ from SBMFT that does not depend on treating 
the nearest neighbour hopping term in the Hamiltonian, 
i^ £ nn , on a different footing than H^ nn . 

Here we instead divide the full Hamiltonian into its 
Heisenberg piece, Hj H1 and its hopping terms, H™ n + 
H? nn . We first focus on treating H Jh in SBMFT. Using 
standard SBMFT, we find that the spinon propagator is 
now given by 



£o differs from £ in that t(x) = 3/8g s (x)JHX instead of 
t(x) — gt(x)to + 3/8g s (x)JHX where to is the bare nearest 
neighbour hopping strength, Jh is the bare Heisenberg 
coupling, g s (x) is the amount Jh is renormalized in the 
Gutzwiller projection, and x — 0.338. We note that at 
zero doping (x = 0), this propagator coincides with the 
SB propagator in Eqn. (2). 

To treat the remaining hopping terms, H™ n +H™ nn , we 
now proceed as we did with H^ nn alone in the main text. 
Using an RPA approximation, we find for this version of 
the YRZ propagator 
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where now 



where 
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^o(k) = — 2t(x)(cos k x + cos k y ). 



(7) e(k) = -2g t (x)to + eM- 

Most importantly, this YRZ propagator retains a line 
of Luttinger zeros along the magnetic BZ or Umklapp 
surface. 



